Genomic epidemiology of SARS-CoV-2 within households in coastal Kenya: a case ascertained cohort study

BackgroundAnalysis of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) genomic sequence data from household infections should aid its detailed epidemiological understanding. Using viral genomic sequence data, we investigated household SARS-CoV-2 transmission and evolution in coastal Kenya households. MethodsWe conducted a case-ascertained cohort study between December 2020 and February 2022 whereby 573 members of 158 households were prospectively monitored for SARS-CoV-2 infection. Households were invited to participate if a member tested SARS-CoV-2 positive or was a contact of a confirmed case. Follow-up visits collected a nasopharyngeal/oropharyngeal (NP/OP) swab on days 1, 4 and 7 for RT-PCR diagnosis. If any of these were positive, further swabs were collected on days 10, 14, 21 and 28. Positive samples with an RT-PCR cycle threshold of <33.0 were subjected to whole genome sequencing followed by phylogenetic analysis. Ancestral state reconstruction was used to determine if multiple viruses had entered households. ResultsOf 2,091 NP/OP swabs that were collected, 375 (17.9%) tested SARS-CoV-2 positive. Viral genome sequences (>80% coverage) were obtained from 208 (55%) positive samples obtained from 61 study households. These genomes fell within 11 Pango lineages and four variants of concern (Alpha, Beta, Delta and Omicron). We estimated 163 putative transmission events involving members of the sequenced households, 40 (25%) of which were intra-household transmission events while 123 (75%) were infections that likely occurred outside the households. Multiple virus introductions (up-to-5) were observed in 28 (47%) households with the 1-month follow-up period. ConclusionsWe show that a considerable proportion of SARS-CoV-2 infections in coastal Kenya occurred outside the household setting. Multiple virus introductions frequently occurred into households within the same infection wave in contrast to observations from high income settings, where single introduction appears to be the norm. Our findings suggests that control of SARS-CoV-2 transmission by household member isolation may be impractical in this setting.

Impact and mitigation of sampling bias to determine viral spread: evaluating discrete phylogeography through CTMC modeling and structured coalescent model approximations

Bayesian phylogeographic inference is a powerful tool in molecular epidemiological studies that enables reconstructing the origin and subsequent geographic spread of pathogens. Such inference is, however, potentially affected by geographic sampling bias. Here, we investigated the impact of sampling bias on the spatiotemporal reconstruction of viral epidemics using Bayesian discrete phylogeographic models and explored different operational strategies to mitigate this impact. We considered the continuous-time Markov chain (CTMC) model and two structured coalescent approximations (BASTA and MASCOT). For each approach, we compared the estimated and simulated spatiotemporal histories in biased and unbiased conditions based on simulated epidemics of rabies virus (RABV) in dogs in Morocco. While the reconstructed spatiotemporal histories were impacted by sampling bias for the three approaches, BASTA and MASCOT reconstructions were also biased when employing unbiased samples. Increasing the number of analyzed genomes led to more robust estimates at low sampling bias for CTMC. Alternative sampling strategies that maximize the spatiotemporal coverage greatly improved the inference at intermediate sampling bias for CTMC, and to a lesser extent, for BASTA and MASCOT. In contrast, allowing for time-varying population sizes in MASCOT resulted in robust inference. We further applied these approaches to two empirical datasets: a RABV dataset from the Philippines and a SARS-CoV-2 dataset describing its early spread across the world. In conclusion, sampling biases are ubiquitous in phylogeographic analyses but may be accommodated by increasing sample size, balancing spatial and temporal composition in the samples, and informing structured coalescent models with reliable case count data.

Comparing the evolutionary dynamics of predominant SARS-CoV-2 virus lineages co-circulating in Mexico

Up to November 2021, over 200 different SARS-CoV-2 lineages circulated in Mexico. To investigate lineage replacement dynamics, we applied a phylodynamic approach to explore the evolutionary trajectories of five dominant lineages that circulated during the first year of the local epidemic. For most lineages, peaks in sampling frequencies coincided with different epidemiological waves of infection in the country. Lineages B.1.1.222 and B.1.1.519 showed comparable dynamics, represented by clades likely originating in Mexico and persisting for over a year. Lineages B.1.1.7, P.1 and B.1.617.2 also displayed similar dynamics, characterized by multiple introduction events leading to a few successful extended local transmission chains that persisted for several months. We further explored viral movements across the country, applied within the largest clades identified (belonging to lineage B.1.617.2). Many clades were located within the south region of the country, suggesting that this area played a key role in the spread of SARS-CoV-2 in Mexico.

SARS-CoV-2 genomic surveillance in Rwanda: Introductions and local transmission of the B.1.617.2 (Delta) variant of concern

The emergence of the SARS-CoV-2 Delta variant of concern (lineage B.1.617.2) in late 2020 resulted in a new wave of infections in many countries across the world, where it often became the dominant lineage in a relatively short amount of time. We here report on a novel genomic surveillance effort in Rwanda in the time period from June to September 2021, leading to 201 SARS-CoV-2 genomes being generated, the majority of which were identified as the Delta variant of concern. We show that in Rwanda, the Delta variant almost completely replaced the previously dominant A.23.1 and B.1.351 (Beta) lineages in a matter of weeks, and led to a tripling of the total number of COVID-19 infections and COVID-19-related fatalities over the course of only three months. We estimate that Delta in Rwanda had an average growth rate advantage of 0.034 (95% CI 0.025-0.045) per day over A.23.1, and of 0.022 (95% CI 0.012-0.032) over B.1.351. Phylogenetic analysis reveals the presence of at least seven local Delta transmission clusters, with two of these clusters occurring close to the border with the Democratic Republic of the Congo, and another cluster close to the border with Tanzania. A smaller Delta cluster of infections also appeared close to the border with Uganda, illustrating the importance of monitoring cross-border traffic to limit the spread between Rwanda and its neighboring countries. We discuss our findings against a background of increased vaccination efforts in Rwanda, and also discuss a number of breakthrough infections identified during our study. Concluding, our study has added an important collection of data to the available genomes for the Eastern Africa region, with the number of Delta infections close to the border with neighboring countries highlighting the need to further strengthen genomic surveillance in the region to obtain a better understanding of the impact of border crossings on lowering the epidemic curve in Rwanda.

Considerable escape of SARS-CoV-2 variant Omicron to antibody neutralization

The SARS-CoV-2 Omicron variant was first identified in November 2021 in Botswana and South Africa1,2. It has in the meantime spread to many countries and is expected to rapidly become dominant worldwide. The lineage is characterized by the presence of about 32 mutations in the Spike, located mostly in the N-terminal domain (NTD) and the receptor binding domain (RBD), which may enhance viral fitness and allow antibody evasion. Here, we isolated an infectious Omicron virus in Belgium, from a traveller returning from Egypt. We examined its sensitivity to 9 monoclonal antibodies (mAbs) clinically approved or in development3, and to antibodies present in 90 sera from COVID-19 vaccine recipients or convalescent individuals. Omicron was totally or partially resistant to neutralization by all mAbs tested. Sera from Pfizer or AstraZeneca vaccine recipients, sampled 5 months after complete vaccination, barely inhibited Omicron. Sera from COVID-19 convalescent patients collected 6 or 12 months post symptoms displayed low or no neutralizing activity against Omicron. Administration of a booster Pfizer dose as well as vaccination of previously infected individuals generated an anti-Omicron neutralizing response, with titers 5 to 31 fold lower against Omicron than against Delta. Thus, Omicron escapes most therapeutic monoclonal antibodies and to a large extent vaccine-elicited antibodies.

Context-specific emergence and growth of the SARS-CoV-2 Delta variant

The Delta variant of concern of SARS-CoV-2 has spread globally causing large outbreaks and resurgences of COVID-19 cases1-3. The emergence of Delta in the UK occurred on the background of a heterogeneous landscape of immunity and relaxation of non-pharmaceutical interventions4,5. Here we analyse 52,992 Delta genomes from England in combination with 93,649 global genomes to reconstruct the emergence of Delta, and quantify its introduction to and regional dissemination across England, in the context of changing travel and social restrictions. Through analysis of human movement, contact tracing, and virus genomic data, we find that the focus of geographic expansion of Delta shifted from India to a more global pattern in early May 2021. In England, Delta lineages were introduced >1,000 times and spread nationally as non-pharmaceutical interventions were relaxed. We find that hotel quarantine for travellers from India reduced onward transmission from importations; however the transmission chains that later dominated the Delta wave in England had been already seeded before restrictions were introduced. In England, increasing inter-regional travel drove Deltas nationwide dissemination, with some cities receiving >2,000 observable lineage introductions from other regions. Subsequently, increased levels of local population mixing, not the number of importations, was associated with faster relative growth of Delta. Among US states, we find that regions that previously experienced large waves also had faster Delta growth rates, and a model including interactions between immunity and human behaviour could accurately predict the rise of Delta there. Deltas invasion dynamics depended on fine scale spatial heterogeneity in immunity and contact patterns and our findings will inform optimal spatial interventions to reduce transmission of current and future VOCs such as Omicron.

Leveraging of SARS-CoV-2 PCR cycle thresholds values (Ct) to forecast COVID-19 trends

IntroductionWe assessed the usefulness of SARS-CoV-2 RT-PCR cycle thresholds (Ct) values trends produced by the LHUB-ULB (a consolidated microbiology laboratory located in Brussels, Belgium) for monitoring the epidemics dynamics at local and national levels and for improving forecasting models. MethodsSARS-CoV-2 RT-PCR Ct values produced from April 1, 2020, to May 15, 2021, were compared with national COVID-19 confirmed cases notifications according to their geographical and time distribution. These Ct values were evaluated against both a phase diagram predicting the number of COVID-19 patients requiring intensive care and an age-structured model estimating COVID-19 prevalence in Belgium. ResultsOver 155,811 RT-PCR performed, 12,799 were positive and 7,910 Ct values were available for analysis. The 14-day median Ct values were negatively correlated with the 14-day mean daily positive tests with a lag of 17 days. In addition, the 14-day mean daily positive tests in LHUB-ULB were strongly correlated with the 14-day mean confirmed cases in the Brussels-Capital and in Belgium with coinciding start, peak and end of the different waves of the epidemic. Ct values decreased concurrently with the forecasted phase-shifts of the diagram. Similarly, the evolution of 14-day median Ct values was negatively correlated with daily estimated prevalence for all age-classes. ConclusionWe provide preliminary evidence that trends of Ct values can help to both follow and predict the epidemics trajectory at local and national levels, underlining that consolidated microbiology laboratories can act as epidemic sensors as they gather data that are representative of the geographical area they serve.

Genomic Sequencing of SARS-COV-2 in Rwanda: evolution and regional dynamics

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), responsible for coronavirus disease 19 (COVID-19), is a single-stranded positive-sense ribonucleic acid (RNA) virus that typically undergoes one to two single nucleotide mutations per month. COVID-19 continues to spread globally, with case fatality and test positivity rates often linked to locally circulating strains of SARS-CoV-2. Furthermore, mutations in this virus, in particular those occurring in the spike protein (involved in the virus binding to the host epithelial cells) have potential implications in current vaccination efforts. In Rwanda, more than twenty thousand cases have been confirmed as of March 14th 2021, with a case fatality rate of 1.4% and test positivity rate of 2.3% while the recovery rate is at 91.9%. Rwanda started its genomic surveillance efforts, taking advantage of pre-existing research projects and partnerships, to ensure early detection of SARS-CoV-2 variants and to potentially contain the spread of variants of concern (VOC). As a result of this initiative, we here present 203 SARS-CoV-2 whole genome sequences analyzed from strains circulating in the country from May 2020 to February 2021. In particular, we report a shift in variant distribution towards the newly emerging sub-lineage A.23.1 that is currently dominating. Furthermore, we report the detection of the first Rwandan cases of the VOCs, B.1.1.7 and B.1.351, among incoming travelers tested at Kigali International Airport. We also discuss the potential impact of COVID-19 control measures established in the country to control the spread of the virus. To assess the importance of viral introductions from neighboring countries and local transmission, we exploit available individual travel history metadata to inform spatio-temporal phylogeographic inference, enabling us to take into account infections from unsampled locations during the time frame of interest. We uncover an important role of neighboring countries in seeding introductions into Rwanda, including those from which no genomic sequences are currently available or that no longer report positive cases. Our results point to the importance of systematically screening all incoming travelers, regardless of the origin of their travels, as well as regional collaborations for durable response to COVID-19.

Early introductions and community transmission of SARS-CoV-2 variant B.1.1.7 in the United States

The emergence and spread of SARS-CoV-2 lineage B.1.1.7, first detected in the United Kingdom, has become a global public health concern because of its increased transmissibility. Over 2500 COVID-19 cases associated with this variant have been detected in the US since December 2020, but the extent of establishment is relatively unknown. Using travel, genomic, and diagnostic data, we highlight the primary ports of entry for B.1.1.7 in the US and locations of possible underreporting of B.1.1.7 cases. Furthermore, we found evidence for many independent B.1.1.7 establishments starting in early December 2020, followed by interstate spread by the end of the month. Finally, we project that B.1.1.7 will be the dominant lineage in many states by mid to late March. Thus, genomic surveillance for B.1.1.7 and other variants urgently needs to be enhanced to better inform the public health response.

Phylogenetic analyses of SARS-CoV-2 B.1.1.7 lineage suggest a single origin followed by multiple exportation events versus convergent evolution

The emergence of new variants of SARS-CoV-2 herald a new phase of the pandemic. This study used state-of-the-art phylodynamic methods to ascertain that the rapid rise of B.1.1.7 "Variant of Concern" most likely occurred by global dispersal rather than convergent evolution from multiple sources.

SARS-CoV-2 genomic characterization and clinical manifestation of the COVID-19 outbreak in Uruguay

COVID-19 is a respiratory illness caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and declared by the World Health Organization a global public health emergency. Among the severe outbreaks across South America, Uruguay has become known for curtailing SARS-CoV-2 exceptionally well. To understand the SARS-CoV-2 introductions, local transmissions, and associations with genomic and clinical parameters in Uruguay, we sequenced the viral genomes of 44 outpatients and inpatients in a private healthcare system in its capital, Montevideo, from March to May 2020. We performed a phylogeographic analysis using sequences from our cohort and other studies that indicate a minimum of 23 independent introductions into Uruguay, resulting in five major transmission clusters. Our data suggest that most introductions resulting in chains of transmission originate from other South American countries, with the earliest seeding of the virus in late February 2020, weeks before the borders were closed to all non-citizens and a partial lockdown implemented. Genetic analyses suggest a dominance of S and G clades (G, GH, GR) that make up >90% of the viral strains in our study. In our cohort, lethal outcome of SARS-CoV-2 infection significantly correlated with arterial hypertension, kidney failure, and ICU admission (FDR < 0.01), but not with any mutation in a structural or non-structural protein, such as the spike D614G mutation. Our study contributes genetic, phylodynamic, and clinical correlation data about the exceptionally well-curbed SARS-CoV-2 outbreak in Uruguay, which furthers the understanding of disease patterns and regional aspects of the pandemic in Latin America.

A world apart: levels and factors of excess mortality due to COVID-19 in care homes. The case of Wallonia - Belgium.

COVID-19 became pandemic in 2020 and causes higher mortality in males (M) than females (F) and among older people. In some countries, like Belgium, more than half of COVID-19 confirmed or suspected deaths occurring in spring 2020 concerned residents of care homes. The high incidence in this population is certainly linked to its peculiar age structure but could also result from its poorer general health condition and/or from a higher contamination through the staff of care homes, while protection equipment and testing capacity were initially limited. To address these issues, we used data from Wallonia (Belgium) to characterize the distribution of death rates among care home institutions, to compare the dynamics of deaths in and outside care homes, and to analyse how age and sex affected COVID-19 death rates inside and outside care homes. We also used annual death rates as a proxy for the health condition of each population. We found that: (1) COVID-19 death rate per institution varied widely from 0{per thousand} to 340{per thousand} (mean 43{per thousand}) and increased both with the size of the institution (number of beds) and with the importance of medical care provided. (2) 65% of COVID-19 deaths in Wallonia concerned residents of care homes where the outbreak started after but at a faster pace than the outbreak seen in the external population. (3) The impact of age on both annual and COVID-19 mortality closely follows exponential laws (i.e. Gompertz law) but mortality was much higher for the population living in care homes where the age effect was lower (mortality rate doubling every 20 years of age increment in care homes, 6 years outside them). (4) Both within and outside care homes, the ratio of M/F death rates was 1.6 for annual mortality but reached 2.0 for COVID-19 mortality, a ratio consistent among both confirmed and suspected COVID-19 deaths. (5) When reported to the annual death rate per sex and age, the COVID-19 relative mortality was little affected by age and reached 24% (M) and 18% (F) of their respective annual rate in nursing homes, while these percentages reduced to 10% (M) and 9% (F) in homes for elderly people (with less medical assistance), and to 5% (M) and 4% (F) outside of care homes. In conclusion, a c. 130x higher COVID-19 mortality rate found in care homes compared to the outside population can be attributed to the near multiplicative combination of: (1) a 11x higher mortality due to the old age of its residents, (2) a 3.8x higher mortality due to the low average health condition of its residents, and (3) probably a 3.5x higher infection rate (1.6x in homes for elderly people) due to the transmission by its staff, a problem more acute in large institutions. Our results highlight that nursing home residents should be treated as a very specific population, both for epidemiological studies and to take preventive measures, due to their extreme vulnerability to COVID-19.

Evolution and epidemic spread of SARS-CoV-2 in Brazil

Brazil currently has one of the fastest growing SARS-CoV-2 epidemics in the world. Due to limited available data, assessments of the impact of non-pharmaceutical interventions (NPIs) on virus transmission and epidemic spread remain challenging. We investigate the impact of NPIs in Brazil using epidemiological, mobility and genomic data. Mobility-driven transmission models for Sao Paulo and Rio de Janeiro cities show that the reproduction number (Rt) reached below 1 following NPIs but slowly increased to values between 1 to 1.3 (1.0-1.6). Genome sequencing of 427 new genomes and analysis of a geographically representative genomic dataset from 21 of the 27 Brazilian states identified >100 international introductions of SARS-CoV-2 in Brazil. We estimate that three clades introduced from Europe emerged between 22 and 27 February 2020, and were already well-established before the implementation of NPIs and travel bans. During this first phase of the epidemic establishment of SARS-CoV-2 in Brazil, we find that the virus spread mostly locally and within-state borders. Despite sharp decreases in national air travel during this period, we detected a 25% increase in the average distance travelled by air passengers during this time period. This coincided with the spread of SARS-CoV-2 from large urban centers to the rest of the country. In conclusion, our results shed light on the role of large and highly connected populated centres in the rapid ignition and establishment of SARS-CoV-2, and provide evidence that current interventions remain insufficient to keep virus transmission under control in Brazil. One Sentence SummaryJoint analysis of genomic, mobility and epidemiological novel data provide unique insight into the spread and transmission of the rapidly evolving epidemic of SARS-CoV-2 in Brazil.

A phylodynamic workflow to rapidly gain insights into the dispersal history and dynamics of SARS-CoV-2 lineages

Since the start of the COVID-19 pandemic, an unprecedented number of genomic sequences of the causative virus (SARS-CoV-2) have been generated and shared with the scientific community. The unparalleled volume of available genetic data presents a unique opportunity to gain real-time insights into the virus transmission during the pandemic, but also a daunting computational hurdle if analysed with gold-standard phylogeographic approaches. We here describe and apply an analytical pipeline that is a compromise between fast and rigorous analytical steps. As a proof of concept, we focus on the Belgium epidemic, with one of the highest spatial density of available SARS-CoV-2 genomes. At the global scale, our analyses confirm the importance of external introduction events in establishing multiple transmission chains in the country. At the country scale, our spatially-explicit phylogeographic analyses highlight that the national lockdown had a relatively low impact on both the lineage dispersal velocity and the long-distance dispersal events within Belgium. Our pipeline has the potential to be quickly applied to other countries or regions, with key benefits in complementing epidemiological analyses in assessing the impact of intervention measures or their progressive easement.